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ABSTRACT 

We study the effect of large-scale magnetic fields on the non-axisymmetric inertial- 
acoustic modes (also called p-modes) trapped in the innermost regions of accretion 
discs around black holes (BHs). These global modes could provide an explanation 
for the high-frequency quasi-periodic oscillations (HFQPOs) observed in BH X-ray 
binaries. There may be observational evidence for the presence of such large-scale 
magnetic fields in the disks since episodic jets are observed in the same spectral state 
when HFQPOs are detected. We find that a large-scale poloidal magnetic field can 
enhance the corotational instability and increase the growth rate of the purely hy¬ 
drodynamic overstable p-modes. In addition, we show that the frequencies of these 
overstable p-modes could be further reduced by such magnetic fields, making them 
agree better with observations. 

Key words: accretion, accretion discs - hydrodynamics - waves - black hole - magnetic 
field 


1 INTRODUCTION 

Black-hole (BH) X-ray binaries exhibit a wide range of vari¬ 
abilities (e.g. van der Klis 2006; Remillard & McClintock 
2006; Yu & Zhang 2013) Of particular interest is the High- 
Frequency Quasi-Periodic Oscillations (HFQPOs; Remillard 
& McClintock 2006; Belloni et al. 2012; Belloni & Stella 
2014). They have frequencies (several tens to a few hun¬ 
dreds Hz) comparable to the orbital frequency at the In¬ 
nermost Stable Circular Orbit (ISCO) around the BH (with 
mass M ~ lOM©) and thus provide a unique probe to study 
accretion flows near BHs and the effects of strong gravity. 
These HFQPOs are only observed in the intermediate spec¬ 
tral state of BH X-ray binaries, when the system transitions 
between the low/hard state (with the X-ray emission dom¬ 
inated by power-law hard photons) and the thermal state 
(with the emission dominated by thermal disk photons). In¬ 
terestingly, it is during the intermediate state when episodic 
jets are observed from the BH X-ray binaries (e.g. Fender et 
al. 2004). Currently, the physical origin of HFQPOs remains 
elusive, and a number of ideas and models have been sug¬ 
gested or explored, including the orbital motion of hot spots 
in the disc (Stella et al. 1999; Schnittman & Bertschinger 
2004; Wellons et al. 2014), nonlinear resonances (Abramow- 
icz & Kluzniak 2001; Abramowicz et al. 2007), and oscilla- 
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tions of finite accretion tori (Rezzolla et al. 2003; Blaes et al. 
2006). A large class of models are based on BH discoseismol- 
ogy, in which HFQPOs are identihed as global oscillations 
modes of the inner accretion discs (see Kato 2001, Ortega- 
Rodriguez et al. 2008 and Lai et al. 2013 for reviews). 

In this paper, we study the effect of large-scale mag¬ 
netic fields on non-axisymmetric (m > 0) p-modes [also 
called inertial-acoustic modes; see Kato (2001) and Wagoner 
(2008) for reviews] trapped in the innermost region of the 
accretion disc around a BH. Lai & Tsang (2009) and Tsang 
& Lai (2009b) showed that these modes, which consist of 
nearly horizontal oscillations with no vertical structure, can 
grow in amplitude due to wave absorption at the corota¬ 
tion resonance (where the wave pattern speed uj/m matches 
the disc rotation rate fl). This overstability requires that 
the disc vortensity, C, = K;^/(2f2E), where k is the radial 
epicyclic frequency and E is the surface density (see Horak 
& Lai 2013 for the full general relativistic version of vorten¬ 
sity), have a positive gradient at the corotation radius Vc 
(Tsang & Lai 2008; see also Narayan et al. 1987). General 
relativity (GR) plays an important role: For a Newtonian 
disc, with B = K oc and a relatively flat E(r) profile, 

dC,/dr < 0, so corotational wave absorption leads to mode 
damping; with GR, however, k is non-monotonic in the in¬ 
ner disc (e.g., for a Schwarzschild BH, k reaches a maximum 
at r = SiGM/c? and goes to zero at risco = QGMjc?), thus 
p-modes with frequencies such that dC,/dr > 0 at are 
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overstable. Linear calculations based on pseudo-Newtonian 
potential (Lai & Tsang 2009; Tsang & Lai 2009b) and on 
full GR (Horak & Lai 2013) show that the mode frequencies 
approximately agree with the observed HFQPO frequencies 
(with the BH mass and spin as constrained/measured by 
observations), although for rapidly rotating BHs, the linear 
mode frequencies are somewhat too large (Horak & Lai 2013; 
see below). Nonlinear simulations show that these overstable 
modes can grow to the large amplitudes and maintain global 
coherence and well-defined frequencies (Fu & Lai 2013). The 
mode growth and saturation can also be enhanced by turbu¬ 
lent viscosity (Miranda, Horak & Lai 2015). Overall, these 
hydrodynamical studies show that non-axisymmetric disc p- 
modes is a promising candidate to explain HFQPOs in BH 
X-ray binaries, although a robust diagnostics remains out of 
reach due to the complexity of the real systems. 

Magnetic fields are likely present in BH accretion discs. 
They may be created as a result of the nonlinear devel¬ 
opment of the magneto-rotational instability (MRI; Bal- 
bus & Hawley 1998). Large-scale poloidal fields may be ad- 
vected inward with the accretion flow, building up signif¬ 
icant strength in the inner disc (e.g., Lubow et al. 1994; 
Lovelace et al. 2009; Gullet & Ogilvie 2012,2013; Gao & 
Spruit 2013). Such large-scale fields can lead to production 
of jets/outflows from accretion discs through the magneto¬ 
centrifugal mechanism (e.g., Blandford & Payne 1982). Mag¬ 
netic fields threading the BH can also lead to relativistic jets 
from the rotating BHs (McKinney et al. 2012). As noted be¬ 
fore, in BH X-ray binaries, episodic jets are observed in the 
intermediate state, and this is the same state during which 
HFQPOs are detected. This suggests that to properly study 
the BH disc oscillation modes, it is necessary to include 
the effect of large-scale magnetic fields. In our model, the 
disc and corona (coupled by a large-scale poloidal magnetic 
field) oscillate together. We suggest that in the intermedi¬ 
ate state, large-scale magnetic fields are created (perhaps 
episodically; cf. Yuan et al. 2009). This would allow the pro¬ 
duction of episodic jets and the disk oscillations to manifest 
as HFQPOs in hard X-rays as observed (Remillard & Mc- 
Clintock 2006; Belloni et al. 2012). Tagger and collaborators 
(Tagger & Pellat 1999; Varniere & Tagger 2002; Tagger & 
Varniere 2006) have developed a similar picture of disc os¬ 
cillations, which they termed accretion-ejection instablity, 
although they focused on the MHD form of Rossby wave in¬ 
stability in the disc (see Lovelace et al. 1999; Yu & Li 2009; 
Yu & Lai 2013). We note that coherent toroidal magnetic 
fields in the disc tend to suppress the corotational instability 
(Fu & Lai 2011). We will focus on poloidal field configura¬ 
tions in this paper. 

Recent study by Horak & Lai (2013) provides a full 
GR corotation instability criterion for disc p-modes. Their 
full GR results of the mode frequencies are qualitatively in 
agreement with the pseudo-Newtonian results (Lai & Tsang 
2009), but indicate that the theoretical frequencies are too 
high compared to observations. The discrepancy is most se¬ 
vere for rapidly spinning BHs (such as GRS 1905-1-105). We 
show in this paper that including large-scale poloidal fields 
in the disc may resolve this discrepency. 

In our model, we consider fluid perturbations that 
have no vertical structure inside the disc (i.e., the vertical 
wavenumber kz = 0). Thus we do not include MRI, which 
generally involves perturbations with finite kz- Since the 


growth of the p-modes is primarily due to in-disc motion, 
it is adequate to consider 2D disc dynamics. Similar setup 
has previously been considered by various authors in differ¬ 
ent contexts (e.g. Spruit et al. 1995; Tagger & Pallet 1999; 
Lizano et al. 2010). 

Our paper is organized as follows. In §2, the basic disc 
equations with large-scale magnetic fields are derived. In §3, 
we present the results of p-modes in thin magnetized discs. 
In §4, we consider the effects of finite disc thickness, and we 
conclude in §5. 


2 BASIC EQUATIONS 


We consider a geometrically thin disc and adopt cylin¬ 
drical coordinate system {r,(j),z). We use the the pseudo- 
Newtonian potential of Paczynski & Wiita (1980) 


$ = - 


GM 
r — rs ’ 


( 1 ) 


where = 2GMjc^ is the Schwarzschild radius. For this 
potential, the free particle (Keplerian) orbital frequency is 



and the radial epicyclic frequency k is given by 


( 2 ) 


K = (3) 

V r-Vs 

The innermost stable circular orbit (ISCO), defined by hd = 
0, is located at risco = 3rs. The epicyclic frequency reaches 
a peak at rmax = (2 -|- y/3)rs. Note that the GR effect plays 
an essential role in the corotational instability of p—modes. 
In Newtonian theory, k = ft ex so dC^/dr < 0 if the 

surface density E is constant; with GR, dC,ldr > 0 for r < 

rmax- 

When large-scale magnetic fields thread the thin con¬ 
ducting disc, the height-integrated mass continuity equation 
and momentum equation read 


-f Vx • (Eu) — 0, (4) 

where Vx is 2D operator acting on the disc plane, E, 
P and u are the surface density, height-integrated pres¬ 
sure and height-averaged velocity, respectively, g = —g{r)r, 
with g — GM/{r — rs)^, is the gravitational acceleration, 
[B]); = B (2 = H) — B( 2 ; = —H) (with H the half-thickness 
of the disc), and we have used = 0. The effect of finite 

disc thickness will be considered in Section 4. 

The unperturbed velocity of an equilibrium disc is 
u = (0,rH,0) (in cylindrical coordinates), with the angu¬ 
lar velocity given by 


- hV = - ^ 


_ ^ 

E dr dr 27rE 


(6) 


where = Br{z = H) — —B^ . The unperturbed disc has 
B^ = 0. Note that Br is nonzero only outside the discs, 
which has different signs above and below the disc. Inside 
the disk, Br is zero (so that the differential rotation of the 
disc would not lead to generation of B,f, inside the disc). 
It should be mentioned that the field configuration 


© 2015 RAS, MNRAS 000, [T]-?? 







Inertial-acoustic Oscillations of Black-Hole Accretion Discs 3 


adopted in this paper is highly idealized. A magneto- 
centrifugal disc wind or jet, if present, will certainly 
involve in the background state. For simplicity, 
we assume that the magnetic field outside the disc to 
be a potential held. The disc surface density is assumed 
to have a power law form E oc r~^, where p is the density 
index. We also try various behavior of the surface 
density and the magnetic field, in which the profile 
of f?z/S can be increasing, decreasing, or constant. 
We find that the results are qualitatively similar. 
Without loss of generality, we fix p = 1 throughout 
this paper. 

The dynamics of the disc is coupled with the large scale 
poloidal magnetic field outside the disc, i.e., the disc magne¬ 
tosphere. We now consider small-amplitude perturbations of 
the disc. We assume that all perturbed quantities have the 
dependence exp{im(f> — iojt), where m = 1, 2, ■ • ■ is the az¬ 
imuthal wave number and uj is the complex frequency. Then 
the linearized fluid perturbation equations become 


—ioj (5S = —Vj_ • (E (5u), 

—iuj Sur — 2n Su^ = 5h - 5B^ 

or 27rE 




Su, + ^6ur = -'-^Sh+^ 5B +, 


(7) 

( 8 ) 
(9) 


where uj = lo — mQ. is the wave frequency in the rotating 
frame of the fluid, and 5h = 5P/T, is the enthalpy pertur¬ 
bation (we assume barotropic discs). The magnetic field in¬ 
duction equation in the disc reads 


— ioj 5Bz = —Vi • (Bz (5u). (10) 

Combining Eqs. 0 and (HQi, we have 



where = Sur/{—iuj) is the Lagrangian displacement. 
Equation (HU can also be derived from {d/dt){Bz/'S) = 0. 
In terms of and 5h, the magnetic field perturbation is 


3Bz = Di^r + D25h, ( 12 ) 


with 

where Cs is the disc sound speed. 

To determine 5Bf and SBf, we assume that the mag¬ 
netic field outside the disc is a potential field (see Spruit et 
al. 1995; Tagger & Pallet 1999). This is equivalent to as¬ 
suming that the Alfven speed above and below the disc is 
sufficiently high that currents are dissipated rapidly (on the 
disc dynamical timescale). Define the “magnetic potential” 
(5 $m outside the disc via 


5B = —sign(2;)V(53>M. (14) 

Then (5 $m satisfies the Poisson equation (Tagger & Pallet 
1999) 

=-2 5Bz(5(a), (15) 

where S{z) is the Dirac delta function. The integral solution 


of (ITSll is 


S^mO) = j <5Bz(r') 1^^ 57 / 2 ( 0 )] dr', (16) 

where a = r'/r, and the Laplace coefficient is defined by 



cos m<j) 

(a^ + eg + 1 — 2a cos 4>)‘^ 


d<j), 


(17) 


with eo the softening parameter (of order the disc aspect 
ratio H/r). The perturbed radial magnetic field at the upper 
disc surface is 


SBf = 

dr 

(s) 


57 ) 2 ( 0 ) + 0 


d5r/2(o) ' 

da 


dr', 


(18) 


and the azimuthal field is SBf = —{im/r)5^M■ 

In our numerical calculations, we use and 5h as the 
basic variables. The perturbation equations are 


d^r 

dr 


^Sh 

dr 


f 2mQ. 1 c^lnE\ . 

—^ + - + —r- O 

\ ru) r dr 


^2(2,2 


5h + 


mr Bz 

r'^dP' 27rE 




-2 2 . B+ V 2 mD 

UJ - K + -I- —oh 

27rE J ruj 

Bz dS^M 2mQ, Bz 

H-5<I>m. 


27rE dr 


ruj 2tvTj 


(19) 


( 20 ) 


To obtain global trapped modes, wave reflection must 
occur at the inner disc boundary rin. To focus on the effect 
of corotational instability, we adopt a simple inner bound¬ 
ary condition, i.e., the radial velocity perturbation vanishes 
at the inner boundary, Sur = 0. This implies a reflect¬ 
ing inner disc edge, for example due to the pres¬ 
ence of a magnetosphere around the black hole (e.g., 
Bisnovatyi-Kogan &: Ruzmaikin 1974, 1976; Igu- 
menshchev et al. 2003; Rothstein &: Lovelace 2008; 
McKinney et al. 2012; see Tsang & Lai 2009a and 
Fu & Lai 2012 for more detailed treatment of the 
magnetosphere-disc interface oscillations). As we are 
interested in the self-excited modes in the inner region of 
the disc, we implement the radiative outer boundary condi¬ 
tions such that waves propagate away from the disc (e.g., Yu 
& Li 2009). The ordinary equations (I19II and (I20II . together 
with the two boundary conditions at the inner and outer 
disk edge, form an eigenvalue problem. Since the eigenfre- 
quency uj is in general complex, equations (O-® are a 
pair of first-order differential equations with complex coef¬ 
ficients which are functions of r. We solve these equations 
using the relaxation method (Press et al. 1992), replacing the 
ODEs by finite-difference equations on a mesh of points cov¬ 
ering the domain of interest (typically risco < r < 2.5risco). 
According to Eqs. (I12II . (1161) and m, both the magnetic 
potential and its derivative can be expressed in terms of the 
linear combination of and Sh. For numerical convenience, 
we calculate the terms in the square bracket in Eqs. unii 
and GHJ and store them for later use. Note that these terms 
are computed only once and can be used repetitively. The 
wave equations (0 and can be cast in a matrix form 
that only deals with the variables d^r and dh. Standard 
relaxation scheme can be applied to the resulting matrix. 
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Figure 1. Rotation curve for different value of (3 when Bf = 
Bz- With the increase of magnetic stress, the angular velocity is 
reduced. 


We use uniform grid points in our calculations. The grid 
point number is typically chosen to be 350. The relaxation 
method requires an initial trial solution that is improved 
by the Newton—Raphson scheme. After iterations the ini¬ 
tial trial solution converges to the eigenfunction of the the 
two-point boundary eigenvalue problem. 


3 RESULTS FOR THIN DISCS 



10 ' 

P 



Figure 2. The growth rate (in units of Qko = ^isco) (upper 
panel) and the real frequency (lower panel) of the m = 2 p- 
mode as a function of the plasma (3 parameter for thin discs [with 
the corresponding vertical magnetic field strength Bz given by 
Eq. (I 22 I 1 I. The blue dashed line is the case of B^ = 0, and the 
red solid line B^ = Bz. 


For our numerical calculations, we adopt the disc sound 
speed Cs = O.YvQk). The magnitude of Bz is specified by 
the dimensionless ratio 

(Eon2ro)i/2’ (21) 

where the subscript “0” implies that the quantities evaluated 
at r = risco- The corresponding plasma parameter in the 
disc at r = risco is 


Po 


STrpCs 


r = 


^isco 


47rRo 

ro 


R 


-2 
zO ’ 


( 22 ) 


where we have used H = Cs/^k — Cs/fl. The magnetic 
field is chosen such the plasma (3 is equal to /3o at r = 
risco throughout the disc. Note that magnetic field Bz in the 
disc varies approximately as Bz oc where p is the 

surface density power law index. We solve for the equilibrium 
rotation profile using Eq. ((6|. Note that when B^ — 0, 
the rotation profile of the disc is unaffected by Bz- 
When B^ — Bz, the equilibrium rotation profile of 
the disc is changed from the nonmagnetic disc. In 
Fig. [l] we show the rotation curve in the case of 
B^ = Bz for different value of p. We clearly see the 
reduction in the angular velocity of the equilibrium 
disc due to the outward BzB^ stress on the disc. 

The upper panel of Figure 2 shows the growth rate of 
the p-modes for different magnetic field strengths. When 
B+ = 0, we see that with the increase of magnetic field, 
the p-mode growth rate varies in a non-monotonic way. The 
growth rate first increases and then decreases. At /? ~ 0.5, 



r/r. 

I SCO 


Figure 4. The angular momentum flux associated with the p- 
mode in disc with B^ = 0. The blue dotted line is for the weakly 
magnetized disc (/3 600), and the red dashed line and black 

solid line are for more strongly magnetized disc {0 = 3.4 and 
0.6). In all cases, the eigenfunction is normalized so that the 
maximum velocity perturbation, 6ur, equals unity. The higher 
angular momentum flux for the more strongly magnetized case is 
consistent with the larger mode growth rate. The vertical lines 
denote the corotation resonances, which are located at 
Tc = 1.19,1.21,1.26, respectively. 
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r/r. r/r. 

ISCO ISCO 

Figure 3. Eigenfunctions of the m = 2 the p-modes trapped in the inner most region in the disc. The left panels show the weakly 
magnetized case (/3 ~ 600); the difference with the nonmagnetic case is essentially indistinguishable. The right panels have /? ~ 0.6 
(strongly magnetized). The top panels show 5ur (the dotted line is for the real part, the dashed line the imaginary part, and the solid 
line the absolute value), and the bottom panels show 5u^. The corotation resonances are marked by the vertical lines, which 
are located at = 1.19,1.26, respectively. 


the growth rate is enhanced by a factor of 3 compared to 
the — 0 value. When Bf = B^, the overall behav¬ 
ior of the growth rate is similar to the case of Bf = 0, 
but the maximum growth rate is about 30% higher than 
disc with Bf = 0. The mode frequencies are shown in the 
lower panel in Fig. 2. For the disc with Bf — 0, the mode 
frequency decreases monotonically from ojr — 0.7m^lKO to 
iJr — O.QmVlKO with the increase of magnetic field. For disc 
with Bf = Bz, the frequency can be reduced by a factor 
about 2, changing from ojr — O.Tmriifo to cjJt — 0.35mflK0- 

Figure [3] gives some examples of the m = 2 eigenfunc¬ 
tions of p-modes trapped in the innermost region of the 
disc. The left panels show the case with very weak magnetic 
fields, and the the right panels show the magnetized case 
with Bf = Bz = 0.78. The amplitudes of the eigenfunc¬ 
tions are normalized so that the maximum absolute value 
of the radial velocity perturbation I^Ur-l equals unity (this 
maximum occurs at r ~ ro). 

To understand the origin of the enhanced p-mode 
growth rate for magnetized discs, we show in Fig. U the 
angular momentum flux associated with the eigenmode for 
the weakly magnetized disc (plasma /3 ~ 600) and for the 
more strongly magnetized discs (plasma P = 3.4 and 0.6). 
The angular momentum flux F{r) across a cylinder of radius 
r is given by (e.g., Goldreich & Tremaine 1979) 

F{r) = J T,5ur5uti,d(j)^ = yrEr^K (&Ur5uf) (23) 

where () designates time average and the superscript * de¬ 
notes complex conjugate. Note that waves carry negative 
(positive) angular momentum inside (outside) the corota¬ 


tion. The net positive (outward) angular momentum flux 
F{r) around the corotation indicates the growth of the p- 
modes. Higher angular momentum fluxes imply higher in¬ 
stability growth rates (see Fig. [4]). 


4 EFFECTS OF FINITE DISC THICKNESS 
4.1 Model Equations 


When finite thickness of the disc is considered, an “internal” 
magnetic force term should be added to the right-hand side 
of Eq. ((5]): 


f 




+ -(B.Vx)B 


(24) 


Obviously, to include the 3D effect rigorously would require 
examining the vertical stratification of the density and mag¬ 
netic field inside the disc - this is beyond the scope of this 
paper. Here we consider a simple model where the internal 
density of the disc is assumed to be independent of so 
that 


S = 2pH, P = 2pH. 

Then the internal magnetic force simplifies to 


f = 


47rp 


-Vj 


B^ 


-f- (B- Vu)B 


(25) 


(26) 


We assume that only vertical magnetic field exists inside the 
unperturbed disc. The equilibrium rotational profile is then 
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determined by 

= . (27) 

E dr 27rE 27rE \ dr J 

To derive the modified perturbation equations including (5f, 
it is convenient to define a new perturbation variable 511 in 
place of 5h-. 



^rr — I BzSBz 2 BzSBz f ^ 

+ (28) 

After some algebra, the final disc perturbation equations can 
be written in the following form: 

d£r (2mQ, 1 d In E „ \ ^ 

- I —^ + - + — +Di]^r 


dr 


\ ruj 


+ 


d 

dr 


m 

r'^uj'^ 

2 2 


dr 


-1>5 5n + 


m Bz 
r'^dS^ 27rE 


5$a 


(29) 


Bt 


5n = + i:>6 e 


+ 


27rE 

(2mQ ^ dlnn 

-r- + Dr- 
V ruj dr 


jn 


p 



p 


Figure 5. Growth rate (in units of f^KO = flisco) and frequency 
of the unstable m = 2 p-mode as a function of the plasma p 
parameter for discs with finite thickness. The blue dashed line is 
for the case with = 0, while the red solid line for the case of 


2mf2 Bz Bz dS^M 

rib 27rE ^ 27rE dr 


(30) 


ST =Bz 


where = Di, and 


D4 


- d( 
ci + ddr y^BzJ' 


Dr. 


De 


Dr 


1 

c^+cl’ 

Da, 

Dr. 


1 ‘^(P+ 
p dr 

1 d{P+WdBz) 
p dr 


Arrp’ 


(31) 

(32) 

(33) 

(34) 


4.2 Results 

Finite disc thickness has two-fold effects on the p-mode coro¬ 
tation instability : the first is the change in the equilibrium 
rotation profile and the second is the direct effect of 5f. With 
these two effects included, we calculate the eigenmodes of 
the disc using the same method as in Section 3. The re¬ 
sults are shown in Figure [5] The red solid line shows the 
case with B^ = Bz and the blue dashed line shows the case 
with B^ = 0. We find that for the parameter considered 
(Cs = Q.lrQ,K, or Hjr = 0.1), the p-mode growth rate is 
decreased compared to the Hjr —>■ 0 limit. More specifi¬ 
cally, the maximum growth rate is about 25% — 30% lower. 
However for sufficiently strongly-magnetized disc (plasma 
P ~ 0.4), especially for disc with radial magnetic field in the 
magnetosphere, the oscillation frequency of the overstable 
mode can be reduced further, changing from Wr — O.TmflKO 
to cjr — 0.2mflK0, which is reduced by about of factor 3 
compared to the pure hydrodynamic mode. 

As noted in Section 1, full GR calculations of hydrody¬ 
namic disc p-modes show that, for accretion discs around ex¬ 
treme Kerr BHs, the oscillation frequency is close to mQ,K 0 
(Horak & Lai 2013). However, observations reveal that the 
frequencies of HFQPOs are usually smaller than this value 
by a factor of 2-3. This discrepancy indicates that pure hy¬ 
drodynamic disc p-modes without magnetic field can not 


explain HFQPOs in X-ray binaries, particularly for high- 
spin sources (such as GRS 1905-1-105). When the large-scale 
poloidal magnetic field is included, this discrepancy can be 
resolved since the reduction in oscillation frequency makes 
the theoretical p-mode prediction more consistent with the 
observed values. 

In addition, we find that the first effect, i.e., the change 
of the equilibrium disc rotation profile induced by the 
disk finite thickness, plays a minor role in the disc p-mode 
instability. To make this point more clearly, we artificially 
discard the last term on the right hand side of Equation 
123. We find that results are very close to those with this 
term included (the detailed results are not shown in Fig. 5, 
since they are very close to the curves in Fig. 5). 


5 DISCUSSION AND CONCLUSION 

In this paper we have carried out linear analysis of p-modes 
in BH accretion discs threaded by large-scale magnetic fields. 
These modes reside in the innermost region in the disc, can 
become overstable driven by corotational instability, and 
may be responsible for the observed high-frequency quasi- 
periodic oscillations (HFQPOs) in BH X-ray binaries. Our 
results show that the large-scale magnetic field can increase 
the mode growth rates significantly (by a factor of a few 
for plasma ~ 1 in the disc). It can also reduce the mode 
frequencies, bringing the theoretical values into agreement 
with the observed HFQPO frequencies. 

It is important to note that since episodic jets are ob¬ 
served in the intermediate state of BH X-ray binaries when 
HFQPOs are detected, a proper treatment of large-scale 
magnetic fields is crucial for understanding the disc-jet- 
QPO connections in accreting BH systems (e.g. Fender et 
al. 2004). Although our treatment of the disc magnetic fields 
is somewhat idealized (e.g., toroidal fields, vertical stratifica¬ 
tion and MRI turbulence are not included), it demonstrates 
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the importance of including magnetic fields in studying BH 
disc oscillations and confronting with observations. 

Overall, our study shows that global oscillations of BH 
accretion discs, combining the effects of general relativity 
and magnetic fields, is a promising candidate for under¬ 
standing HFQPOs observed in BH X-ray binaries and simi¬ 
lar oscillations potentially observed in intermediate and su- 
permassive BHs (e.g., Pasham, Strohmayer & Mushotzky 
2014). Important caveats and uncertainties (see Lai et 
al. 2013 for a concise review of various issues) remain in 
our treatment of the physical condition of the innermost ac¬ 
cretion flows around BHs (e.g., the BH magnetosphere - disc 
interface; see Tsang & Lai 2009a; Fu & Lai 2012; Miranda 
et al. 2015). It is hoped that future observations with more 
sensitive X-ray timing telescopes (such as LOFT, see Feroci 
et al. 2014) will shed more light on the origin of variabilities 
of accreting BHs (e.g. Belloni & Stella 2014). 
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